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tained  at  incipient  breakdown  at  the  most  critical  location  is  pro¬ 
portional  to  the  fluid  density  at  that  location.  It  is  shown  that 
as  a  consequence  of  this  fact,  the  electrical  output  can  be  maxi¬ 
mized  by  designing  the  primary  jet  for  an  exit  Mach  number  of  0.71. 

Estimates  are  made  of  the  pump  work  required,  of  mixing  losses  in 
the  ejector  and  of  friction  and  secondary  flow  losses.  The  mathe¬ 
matical  analysis  is  reduced  to  a  fully  non-dimensional  form  and 
the  key  dimensionless  parameters  that  govern  performance  are 
clearly  identified.  A  preliminary  estimate  is  made  of  the  numerical 
values  of  these  parameters  and  the  overall  performance  of  the  system 
is  estimated  on  this  basis 


Unfortunately,  the  results  so  obtained  are  very  pessimistic.  They 
indicate  that  even  at  100  atmospheres  pressure,  electrical  break¬ 
down  so  severely  limits  power  output,  that  it  is  probably  insuffi¬ 
cient  to  cover  the  demand  for  pump  power  and  for  power  to  overcome 
the  various  losses  involved,  let  alone  provide  any  useful  net  out¬ 
put.  The  only  hopeful  note  is  that  the  analysis  so  clearly  pin¬ 
points  the  problem  that  it  might  in  the  end  also  suggest  the  means 
for  surmounting  it. 

A  modified  design  where  the  electrical  conversion  section  is  placed 
after  the  ejector  is  analyzed  in  the  Appendix.  This  modification 
yields  less  pessimistic  results,  but  still  leaves  doubt  as  to  the 
viability  of  this  type  of  power  generator. 
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ABSTRACT 


-V 

This  report  develops  a  detailed  analysis  of  a  type  of 
electrohydrodynamic  power  generator  which  employs  an  ejector 
and  a  so-called  "fluid  flywheel"  as  essential  components. 

The  medium  is  steam  containing  electrically  charged  water 
droplets. 

The  analysis  takes  into  account  the  experimentally  estab¬ 
lished  fact  that  the  maximum  strength  of  the  electrical  field 
that  can  be  sustained  at  incipient  breakdown  at  the  most  criti¬ 
cal  location  is  proportional  to  the  fluid  density  at  that  loca¬ 
tion.  It  is  shown  that  as  a  consequence  of  this  fact,  the 
electrical  output  can  be  maximized  by  designing  the  primary 
jet  for  an  exit  Mach  number  of  0.71. 

Estimates  are  made  of  the  pump  work  required,  of  mixing 
losses  in  the  ejector  and  of  friction  and  secondary  flow 
losses.  The  mathematical  analysis  is  reduced  to  a  fully  non- 
dimensional  form  and  the  key  dimensionless  parameters  that 
govern  performance  are  clearly  identified.  A  preliminary 
estimate  is  made  of  the  numerical  values  of  these  parameters 
and  the  overall  performance  of  the  system  is  estimated  on  this 
basis.  ^ 

Unfortunately,  the  results  so  obtained  are  very  pessimistic. 
They  indicate  that  even  at  100  atmospheres  pressure,  electrical 
breakdown  so  severely  limits  power  output,  that  it  is  probably 
insufficient  to  cover  the  demand  for  pump  power  and  for  power 
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to  overcome  the  various  losses  involved,  let  alone  provide  any 
useful  net  output.  The  only  hopeful  note  is  that  the  analysis 
so  clearly  pinpoints  the  problem  that  it  might  in  the  end  also 
suggest  the  means  for  surmounting  it. 

A  modified  design  where  the  electrical  conversion  section 
is  placed  after  the  ejector  is  analysed  in  the  Appendix.  This 
modification  yields  less  pessimistic  results,  but  still  leaves 


doubt  as  to  the  viability  of  this  type  of  power  generator. 


1 .  Introduction 


The  purpose  of  this  analytical  report  is  to  develop  a  mathe¬ 
matical  model  of  a  proposed  type  of  electrohydrodynamic  (EHD)  power 
generator  which  embodies  an  ejector  as  an  essential  component.  The 
overall  configuration  is  approximately  toroidal  in  form  in  order 
to  help  conserve  the  angular  momentum  of  the  working  fluid  which 
circulates  through  it.  For  this  reason  the  inventor*  of  this  con¬ 
figuration  has  dubbed  it  a  "fluid  flywheel". 

Fig.  1.1  is  a  schematic  diagram  of  the  device.  Superheated 
steam  at  stagnation  pressure  pQ  and  stagnation  temperature  TQ 
is  supplied  to  a  nozzle  which  discharges  into  the  ejector  mixing 
region  at  station  1.  A  secondary  fluid  enters  the  mixing  region  at 
station  2.  The  resulting  mixture  leaves  the  ejector  at  station  3. 
The  ejector  action  creates  a  pressure  rise  across  the  mixing  sec¬ 
tion  which  is  just  sufficient  to  offset  the  effect  of  the  power  out¬ 
put  and  of  the  various  flow  losses  that  occur  in  the  rest  of  the 
circuit. 

The  mixture  of  liquid  and  vapor  which  leaves  station  3  is  cooled 
during  its  return  circuit  sufficiently  to  condense  and  separate  out 
a  mass  flow  of  liquid  equal  to  the  mass  flow  of  steam  which  passes 
through  the  primary  jet. 

The  liquid  withdrawn  from  the  condenser  is  pumped  back  into  the 
boiler  and  superheater  where  it  is  heated  back  to  the  initial  stag¬ 
nation  state  0.  The  primary  fluid  thus  undergoes  a  Rankine  cycle 
where  the  turbine  and  condenser  have  been  replaced  by  an  EHD 
generator. 

*  See  note  in  Section  9,  References. 
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SCHEMATIC  OF  POWER  GENERATOR 


PUMP 


Conditions  in  the  primary  nozzle  are  controlled  in  such  a 
way  as  to  cause  fine  droplets  of  moisture  to  condense  here  and 
electrodes  are  provided  which  electrically  charge  the  liquid  drop¬ 
lets.  These  electric  charges  are  removed  from  the  flow  by  a  grid 
type  collector  located  at  the  ejector  exit,  station  3.  The  moving 
charges  constitute  an  electric  current  which  sets  up  a  correspond¬ 
ing  electric  field.  The  electric  forces  acting  on  the  charged 
particles  tend  to  oppose  the  motion.  However,  the  size  and  charge 
of  the  particles  is  controlled  in  such  a  way  that  they  have  negli¬ 
gible  mobility  with  respect  to  the  surrounding  fluid.  Thus  the 
moving  gas  does  work  on  the  fluid  particles  in  moving  them  down¬ 
stream  against  the  resistance  of  the  electrical  forces.  In  per¬ 
forming  this  work,  the  gas  stream  undergoes  a  corresponding  de¬ 
crease  in  enthalpy .  The  work  done  by  the  gas  against  the  elec¬ 
trical  forces  sets  up  an  electrical  potential  difference  between 
the  charging  electrode  and  the  collector  grid.  If  these  two 
terminals  be  connected  by  an  external  electrical  circuit,  that 
circuit  will  therefore  deliver  useful  electrical  power.  The 
only  turbo  device  required  is  a  small  pump. 

The  question  that  arises  is  whether  such  a  hypothetical  scheme 
can  be  made  practically  effective.  There  are  various  formidable 
obstacles  to  overcome,  including  the  achievement  of  negligible 
mobility.  It  is  therefore  important  to  have  a  complete  quan¬ 
titative  model  of  all  cycle  components  to  help  guide  research 
and  development  efforts. 
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Such  a  model  has  now  been  developed  for  the  fluid  flywheel 
concept  and  is  presented  in  this  report.  The  model  shows  that  the 
performance  that  can  be  attained  is  governed  by  just  a  few  basic 


dimensionless  parameters.  It  identifies  these  parameters  and 
provides  some  basis  for  estimating  their  probable  magnitudes. 

By  identifying  the  key  parameters  and  showing  their  effect  on 
overall  performance,  this  analysis  permits  current  experimental 
research  efforts  to  be  redirected  along  more  fruitful  lines. 

This  report  is  the  result  of  research  sponsored  by  the  Depart¬ 
ment  of  Energy. 


2.  Basic  Electrical  Performance 

It  is  convenient  in  the  present  analysis  to  introduce  a 
parameter  a  which  is  used  to  distinguish  between  the  two  dis¬ 
tinct  cases  of  positively  and  negatively  charged  particles. 
Specifically,  we  set  a  »  +1  for  the  case  of  positively  charged 
particles,  and  we  set  c  =  -1  for  the  case  of  negatively  charged 
particles.  It  is  then  appropriate  to  denote  the  electrical  charge 
per  unit  mass  by  the  product  cq  where  q  has  the  units  of 
coulombs/kg  and  is  always  positive  by  definition. 

The  electrical  working  section  of  the  generator  coincides 
with  the  mixing  region  of  the  ejector  and,  assuming  polar  sym¬ 
metry,  the  distribution  of  q  is  essentially  two  dimensional 
over  this  region.  However,  our  purpose  in  the  present  section  is 
to  develop  certain  key  dimensionless  parameters  which  depend 
primarily  on  how  q  varies  with  respect  to  the  axial  coordinate 
z.  These  parameters  are  not  sensitive  to  variation  with  respect 
to  the  radial  coordinate  r.  Under  these  circumstances  it  is 
permissible  to  simplify  the  analysis  in  the  following  way:  At 
each  axial  station  z  that  lies  within  the  mixing  region,  we 
subdivide  the  total  cross-sectional  area  A  into  an  inner  circular 
region  of  area  A'  and  an  outer  concentric  annulus  of  area  A" 
such  that 

A'  +  A"  =  A  =  A^  *  constant  along  channel  (2.1) 

While  fluid  density  P  (kg/m)  is  treated  as  essentially  uniform 
over  the  entire  mixing  region,  the  axial  velocity  (m/sec)  is 
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assumed  to  have  the  value  V'  over  the  inner  area  A'  and  the 


value  V"  =*  V2  over  the  outer  annular  area  A".  Similarly,  charge 
density  (coulomb/kg)  is  assumed  to  have  the  value  q'  over  the 
inner  area  A'  and  the  value  q"  =  0  over  the  outer  annular 
area  A".  Moreover,  the  quantities  A',  V'  and  q'  are  regarded 
as  functions  of  coordinate  z  which  vary  monotonically  from  their 
initial  values  A^,  V^,  q^  at  z  =  0  to  their  final  values 
A^,  V^,  q^  at  z  =  1.  The  exact  form  of  these  functions  is  not 
of  immediate  importance  in  the  present  development. 

It  is  useful  to  introduce  the  following  auxiliary  definitions 
at  this  point  on 


l  =  length  of  working  section,  m 

V  =  magnitude  of  voltage  change  across  working  section, 
positive  by  definition,  volts 

e  =  =  permittivity  of  free  space 

=  8.854  x  10“12  farad/m 

$  =  electric  potential,  volts 


f 

? 

g 


^  =  dimensionless  electric  potential 
2 

j  -  dimensionless  axial  coordinate 

q 1 

j*  =  dimensionless  charge/mass  ratio 


Poisson's  equation  and  its  first  two  integrals  can  now  be 
represented  in  the  following  one  dimensional  version: 


V  /dif\ 
vdz2/  l2  \dC2/ 


a 


g(C) 


(2.2) 


2-2 


(2.3) 


0  '  2  ?  C 

(♦-♦.)  -  V(f-f0)  *  (||)  U  -  °P"e-  f  f  9(C)dsdC 


In  these  equations,  the  quantity  ^|  is  the  negative  of  the 
electrical  field  strength  (volts/m) .  This  quantity  attains  its 
maximum  magnitude  at  station  z  =  0,  which  is  the  critical  loca¬ 
tion  at  which  electrical  breakdown  first  occurs. 

It  has  been  established  by  experiment.  Ref.  (1) ,  that  over 
a  certain  range  of  pressures,  the  breakdown  strength  of  a  given 
medium  is  well  approximated  by  the  following  expression 

(ft).-  °[Co  +  cbr°i}  <2-5> 


where  R  ( joules/kg°K)  is  the  gas  constant  of  the  medium  and  where 

2 

the  breakdown  constants  C  (volts/m)  and  C„  (m  °K/ coulomb)  are 

O  D 

characteristic  properties  of  the  medium  which  can  be  measured  ex¬ 
perimentally  . 

The  evaluation  of  Eqs  (2.3)  and  (2.4)  at  the  exit  station 
£  =  1  is  of  particular  importance.  In  this  connection  it  is  conven¬ 
ient  to  introduce  the  following  notation: 

1 

<I>1  ~  J  g(S)d?  (2.6) 

0 
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(2.7) 


1  C 

■// 

0  0 


g(C) dcd? 


Notice  that  i ^  and  <1/ ^  are  dimensionless  coefficients  that 
depend  only  on  the  form  of  the  dimensionless  charge/mass  function 
g(0  . 

The  length  l  of  the  working  section  of  an  ideal  EHD  generator 
having  a  given  current  flow  i  (amps)  should  be  so  chosen  as  to 
maximum  the  potential  difference  V  (volts)  and  hence  the  electri¬ 
cal  power  Pe  (watts)  which  is  generated.  This  means  that  the 
dimensionless  potential  f  should  attain  its  maximum  value  at 
the  exit  station,  or  that 


(if)- 0  at  c  ’ 1 

It  also  follows  that 


(2.8) 


(f  -  f0)  ■  a  at  £  =  1 


(2.9) 


Upon  substituting  Eqs  (2.5)  through  (2.9)  into  Eqs  (2.3)  and 
(2.4)  we  can  readily  extract  a  pair  of  useful  relations  from  the 
result.  In  this  connection  it  is  also  convenient  to  introduce 
the  auxiliary  notation 


\^BRP1/  CB  \Pq/ 


1  + 


1 

rpn 


(2.10) 
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Note  that  a  is  dimensionless.  This  quantity  is  normally  very  small 
compared  with  unity. 

Using  this  notation  and  procedure,  we  may  now  summarize  Eqs 
(2.3)  and  (2.4)  in  the  form 


*1* 


(1+q) 

♦l 


eCgR 


(2.11) 


y1  ”  ^J(l+a)CBRPi£ 


(2.12) 


The  current  flow  i  (amps)  through  the  channel  may  be  written 

i  =  cq1P1V1A1  (2.13) 

The  gross  electrical  power  output  Pg  (watts)  is  now  given 
by  the  simple  expression 

pe  =  [gV]i  (2.14) 

In  order  to  reduce  our  final  result  into  a  more  useful  form 
we  also  write 

v!  =  (2.15) 

where 

V1  =  axial  velocity  to  nozzle  exit,  m/sec 
a1  =  sonic  velocity  at  nozzle  exit,  m/sec 
=  Mach  number  at  nozzle  exit 
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Upon  combining  the  last  five  relations  and  rearranging,  we 


obtain  the  important  result 


Pe  "  lif1  *  J^d^^cl^pjai^ 


(2.16) 


For  the  purpose  of  non-dimensionalizing  this  result,  it  is 
useful  to  introduce  the  following  additional  dimensional  reference 
parameters 


aQ  =  sonic  velocity  corresponding  to  stagnation  conditions 
at  inlet  to  primary  jet,  m/sec 

p  =  density  corresponding  to  the  above  stagnation  conditions, 
°  kg/m 3 


=  throat  area  of  primary  jet,  m^ 


,2-2.2. 


Now  dividing  Eq  (2.16)  through  by  eCgR  pQaQAt 
gives 


^CBRWt, 


=  r  (l  -  2r- )  |(l+a)  2 


vpq /  \ao/  1\Aty 


(2.17) 


Let  us  denote  the  bracketed  quantity  on  the  right  in  the 
following  way 


(Cl+ct)  2  (-1)  Li  M,  Li 
Vo/  Vo/  l\At/ 


>  a  (1+a)  ^i|/3  F  (M1) 


(2.18) 


where  i|>3  is  a  certain  normalizing  constant  which  will  be  defined 
presently. 


In  order  to  analyze  the  function  defined  by  Eq  (2.18),  we  treat 
the  expansion  through  the  primary  nozzle  in  the  usual  way  as  an 
isentropic  process  of  constant  total  enthalpy.  The  medium  is 
treated  as  a  perfect  gas  for  which  y  ,  the  ratio  of  specific  heats , 
is  a  constant. 

The  relations  which  govern  such  a  nozzle  flow  are  well  known 
and  will  therefore  be  used  freely  here  without  detailed  deriva¬ 
tion.  Note  that  there  are  two  distinct  expressions  for  the  func¬ 
tion  depending  on  whether  is  subsonic  or  supersonic. 

For  the  subsonic  case,  the  throat  of  the  nozzle  coincides  with  the 
exit  station  so  that  (A^/Afc)  equals  unity.  In  the  supersonic 
case  (A^/Afc)  becomes  a  function  of  . 

When  these  details  are  properly  worked  out,  the  following 
solution  is  obtained. 


(2.19) 
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By  differentiating  Eq  (2.19),  it  can  be  shown  that  the  function 
passes  through  its  maximum  value  precisely  at 


M 


lopt 


-  i  0.71 

v/5 


(2.21) 


Note,  incidentally,  that  this  result  happens  to  be  independent 
of  y[ 

2 

On  the  other  hand  the  function  (1  +  a)  ^FfM^)  does  not 
necessarily  pass  through  its  maximum  at  precisely  this  same  Mach 
number  because,  as  may  be  seen  from  Eq  (2.10),  the  quantity  a 
is  itself  a  function  of  density  p  at  station  1,  and  hence  of 
M^.  Nevertheless,  a  is  very  small  compared  with  unity,  and  its 
variation  will  therefore  have  an  effect  on  so  small  that 

it  can  safely  be  neglected. 

As  a  normalizing  condition  we  stipulate  that  the  maximum 
value  of  the  function  F(M^)  shall  be  unity.  Thus 


Fmax  -  F<Mlopt>  *  1 


(2.22) 


This  constraint  now  fixes  the  magnitude  of  the  normalizing  con¬ 
stant  which  works  out  to  be 


mk  W 


(Y+3) 

177% 


(2.23) 
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For  example,  for  y  =  1.3  this  formula  gives 


U>3  =  0.4211 


(2.24) 


The  dimensionless  coefficients  <Ji^,  ^  and  i|>3  and  the  factor 
2 

(1  +  a)  can  next  be  combined  into  an  overall  dimensionless  coef¬ 
ficient  according  to  the  relation 


♦4 


(l+o) 


2 


(2.25) 


With  this  notation,  the  electrical  power  output  can  be  expres¬ 
sed  in  the  following  form,  namely. 


— T”3 — 5 — 

■eCB  R  Pc  aoA‘ 


*4f(m1) 


(2.26) 


where  the  function  FfM^)  is  defined  by  Eqs  (2.19),  (2.20) 
and  (2.23)  . 

The  function  F(M^)  is  listed  in  Table  2.1  and  plotted  in 
Fig  2.1.  These  results  are  for  y  =  1.3  which  is  a  typical  value 
often  used  for  steam.  Notice  that  F(M^)  reaches  its  peak  value 
of  unity  at  an  optimum  exit  Mach  number  of  approximately  0.71. 
Observe  that  appreciable  deviations  from  this  optimum  exit  condi¬ 
tion  cause  heavy  performance  losses . 

Under  these  conditions,  the  subsequent  analysis  will  be  re¬ 
stricted  specifically  to  the  optimum  exit  Mach  number  =  0.71 
as  derived  above.  Consequently,  Eq  (2.25)  can  then  be  simplified 


TABLE  2.1 
Y  =  1.3 


0.0 

F  (M^) 

0.0 

0.20 

0.4550 

0.40 

0.8014 

0.60 

0.9774 

0.71 

1.0000 

0.80 

0.9849 

1.00 

0.8722 

1.20 

0.7241 

1.40 

0.5886 

1.60 

0.4704 

1.80 

0.3711 

2 .00 

0.2901 

2.20 

0.2253 

2.40 

0.1744 

2.60 

0.1347 

2.80 

0.1041 

3.00 

0.0805 

3.20 

0.0625 

3.40 

0.0486 

3.60 

0.0380 

3.80 

0.0298 

4.00 

0.0235 

4.20 

0.0186 

4.40 

0.0148 

4.60 

0.0119 

4.80 

0.0095 

5.00 

0.0077 

to  the  form 


Pe 


|eCB-R  p 


“2 - j- 

0  aoAt/ 


(2.27) 


In  order  to  get  the  greatest  benefit  from  Eq  (2.27)  we  should 
know  the  value  of  the  power  coefficient  .  This  coefficient 
is  defined  by  Eqs  (2.6),  (2.7),  (2.23)  and  (2.25).  It  may  be  seen 
from  Eqs  (2.6)  and  (2.7)  that  in  order  to  calculate  the  quantities 
^  and  i|>2  that  ultimately  fix  the  initially  unknown  coefficient 

i|>4,  we  must  know  the  exact  dimensionless  charge/mass  distribution 
as  fixed  by  the  function  g(?)  .  Unfortunately  this  function  re¬ 
mains  unknown  so  that  an  exact  solution  is  not  possible  within  the 
scope  of  the  present  simplified  analysis. 

Nevertheless ,  we  do  have  some  approximate  information  about 
this  function  which  is  adequate  to  provide  a  first  order  estimate 
of  Let  us  call  this  ^4Q  .  In  order  to  establish  such  an 

estimate  or  reference  value  it  is  only  necessary  to  assume  a 
qualitatively  correct  and  analytically  convenient  form  for  this 
function. 

Note  that 

g  =  1  at  ;  =  0  (2.28) 

„  _  _  g3 

g  g3  q^  at  {  =  1  (2.29) 

Moreover  g(£)  varies  monotonically  between  these  two  end 
values  and  it  should  also  satisfy  the  condition,  at  least 
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approximately,  that 


It  is  also  timely  to  note  here  that  the  current  i  evaluated 
at  the  two  stations  1  and  3  gives  the  relation 


qlplAlVl  =  c5-j*33A3V3  (2.31) 

From  this  we  readily  infer  that 

a  _  _/plAlvl\ 

3  ~  ql  \p3A3V3/  ~  X  (2.32) 

where  x  is  the  dimensionless  mass  flow  ratio  of  the  ejector. 

A  simple  analytical  form  for  the  reference  function  g(S) 
that  satisfies  the  foregoing  constraints  is 


g  (?)  =  x+  (1-  x)  (l-£)  2  (2.33) 

Upon  substituting  this  into  Eqs  (2.6)  and  (2.7)  and  inte¬ 
grating  we  obtain  the  expressions 


♦  l  =  j  (1+2X) 

<l>2  =  J  (1+x) 


2-13 


Then  using  Eq  (2.25)  we  obtain  the  desired  reference  value  as 

where  constants  ip3  and  a  are  still  defined  by  Eqs  (2.23)  and 
(2.10)  as  before. 

Of  course  ip4Q  is  only  a  reference,  not  the  true  value  i p4 
itself.  However,  we  may  relate  these  two  quantities  as  follows. 


*4  =  u*40 


(1+a) 


2 


(l+5x) 

(l+2x) 


2 


(2.37) 


Here  u  is  regarded  as  an  empirical  correction  factor  of  order 
unity.  It  is  to  be  expected  that  whereas  \p4  itself  may  tend  to 
vary  over  a  certain  range,  y  should  be  more  nearly  constant. 

We  may  now  combine  Eqs  (2.27)  and  (2.37)  in  the  form 


P 

_ e_ 

2„2 "  2 


'eCB  R  po  aoA* 


y 


(1+a) 


2 


(l+5x) 
(l+2x)  2 


(2.38) 


The  quantities  i|>3,  a  and  x  which  appear  in  Eq  (2.38)  are 
defined  by  Eqs  (2.23),  (2.10)  and  (2.32),  respectively.  The  coef¬ 
ficient  y  must  be  determined  by  experiment  or  by  more  advanced 
analysis.  However,  for  purposes  of  order  of  magnitude  estimates 
it  is  permissible  to  take 


y  =  1 


(2.39) 
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In  analyzing  the  mixing  and  friction  losses  in  the  ejector 
and  the  friction  and  secondary  losses  in  other  parts  of  the  sys¬ 
tem,  it  will  prove  useful  to  employ  as  dimensional  reference 
quantities  the  fluid  density  p=p1®p2=p3  in  the  ejector, 
the  exit  area  A^  from  the  ejector  and  the  average  axial 
velocity  V3  at  this  station.  Thus  the  reference  mass  flow  rate 

becomes  p^A^V^  (kg/sec) ,  the  reference  kinetic  energy  per  unit 
2 

mass  becomes  V^/2  (joule/kg)  and  the  reference  power  becomes 
2 

( p jA^V^)  (V^/2)  (watts) .  In  this  connection  it  is  appropriate  to 
define  a  dimensionless  gross  electrical  power  coefficient  cQ 
in  the  following  way,  namely. 


P 

C  =  - - -  (3.1) 

[P3A3V3]  [v32/2] 

It  is  also  necessary  to  establish  suitable  links  between  the 
above  reference  quantities  and  certain  corresponding  quantities 
at  station  1,  the  exit  from  the  primary  jet.  In  this  connection 
recall  that 


/*l\  /  PlAlVl\ 
\&J  \p3A3V3/ 


dimensionless  mass  flow  ratio 


(3.2) 


We  also  elect  to  introduce  the  following  definition 


quantities,  namely. 


(3.6) 


and 


6 


=  < 


2  2 

Y£cb  R  Pq 


(3.7) 


Recall  that  the  quantities  involved  in  Eq  (3.6)  all  correspond 


lopt 


l/v2. 


Hence  ^  turns  out  to  be  a  constant 


to  the  value  M. 


which  depends  only  on  y  and  a. 
we  find  that 


+  a) 


When  the  details  are  worked  out 


(3.8) 


For  example,  for  y  =1.3  this  gives 


4>  =  1.32(1  +  a)2 


(3.9) 


Since  the  medium  is  treated  as  a  perfect  gas,  the  crucially 
important  parameter  defined  by  Eq  (3.7)  can  be  further  reduced  to 
the  form 


S  = 


(3.10) 


This  result  shows  the  importance  of  maintaining  inlet  stagna¬ 
tion  pressure  pQ  as  high  as  possible  and  inlet  stagnation  tempera¬ 
ture  Tq  as  low  as  possible.  However,  if  condensation  is  used  to 
produce  droplets,  Tq  is  restricted  to  that  range  of  temperatures 
which  produces  optimum  droplet  size.  One  way  around  this  restriction 
is  as  follows:  Let  the  carrier  fluid  be  changed  from  steam  to  air 
(or  to  some  other  gas  which  is  noncondensible  at  ordinary  tempera¬ 
tures)  .  Let  charged  liquid  droplets  or  solid  particles  of  proper 
size  be  injected  into  the  carrier  fluid.  This  permits  Tq  to  be 
sharply  reduced  whereupon  6  is  correspondingly  increased.  Of 
course  this  alternative  may  involve  its  own  characteristic  draw¬ 
backs  as  well,  but  these  will  not  be  discussed  further  at  this 
point. 
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While  the  introduction  of  the  separate  quantities  u,  ip  and  6 
follows  naturally  from  Eq  (3.5),  we  note  that  these  quantities 
always  occur  in  the  form  of  the  product  .  Hence  it  will 

simplify  the  subsequent  development  to  set 


u  =*  ytfiB 

whereupon  Eq  (3.5)  reduces  to 


(3.11) 


c 


e 


U  +  5X) 
1+2X) 2 


(3.12) 


In  the  subsequent  analysis  of  ejector  performance,  it  will  be 
shown  how  the  two  fundamental  parameters  x  and  y  govern  pressure 
rise  through  the  ejector,  mixing  losses  and  so  on.  Eq  (3.12) 
therefore  ties  in  gross  electrical  power  output  with  these  other 
effects  on  the  same  consistent  overall  basis.  In  particular,  it 
is  important  that  the  gross  electrical  power  coefficient  ce  as 


expressed  by  this  equation  should  exceed  the  sum  of  all  the 
losses  in  the  system  by  as  wide  a  margin  as  possible. 
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4.  Ejector  Performance 

Recall  the  notation  introduced  earlier,  namely, 

(4.1) 

*  x  *  mass  flow  ratio 

(4.2) 

=  y  =  velocity  ratio 

As  a  basis  of  comparison  with  the  actual  physical  ejector, 
consider  a  hypothetical  ideal  ejector  which  satisfies  the  following 
constraints.  Firstly,  the  fluid  density  is  uniform  throughout 
the  ejector  so  that  p1=p2=p3=p=  constant.  Secondly,  pres¬ 
sure  is  uniform  across  the  entire  inlet  so  that  p2  =  p1 .  Thirdly, 
the  velocities  are  purely  axial  and  are  uniform  across  sections 
1,  2  and  3.  Fourthly,  wall  friction  is  negligible. 

Imagine  this  ideal  ejector  to  be  in  operation  with  all  elec¬ 
trical  circuits  turned  off.  The  ideal  pressure  (p^*  -  p^)  that 
would  occur  under  these  conditions  may  be  expressed  by  means  of 
the  dimensionless  coefficient 


tg3*  -  Pi} 

w 


ideal  pressure  rise  coefficient 
(for  given  values  of  x  and  y) 


(4.3) 


Of  course  the  real  physical  unit  operating  under  realistic 
conditions  with  electrical  power  output  connected  to  a  load,  will 
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produce  some  smaller  pressure  rise  (p3  -  p^)  at  the  same  values 
of  x  and  y.  This  actual  pressure  rise  may  be  expressed  by  the 
coefficient 


(P3  ~  PL) 

^V32 


actual  pressure  rise  coefficient 
(for  given  values  of  x  and  y) 


(4.4) 


The  difference  (cp*  “  c  )  accounted  for  in  part  by  the  gross 
electrical  power  output  as  expressed  by  the  previously  analyzed 
coefficient  ce  and  in  part  by  additional  mixing  and  friction 

losses  that  characterize  the  real  system. 

The  mass  flow  rates  across  sections  1,  2  and  3  may  be  written 


=  PA1V1  =  xpA3V3 

(4.5) 

=*  pA2V2  =  (1-X)  A3V3 

(4.6) 

*  pA3V3 

(4.7) 

From  these  relations  and  from  Eqs  (4.1)  and  (4.2)  we  find 
the  two  area  ratios 
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These  ratios  satisfy  the  constant  area  condition 


M  /A2\ 

s) +  4 ■ 1 


(4.10) 


We  can  solve  these  last  three  equations  for  the  velocity 


ratio 


2  \  _  (l-x)y 
'  TjFx) 


(4.11) 


The  momentum  equation  can  now  be  applied  to  the  ideal  ejector 
as  follows 

(p3*  -  P]_) A3  =  pA3V32  j  x(^j  +  (1  -  -  1  (4.12) 


We  next  divide  through  by  and  use  Eqs  (4.2)  and  (4.11) 
to  eliminate  the  velocity  ratios  (V^/V^.  and  (V^/V3)  .  After 
simplification,  this  reduces  to  the  important  result 


i2c(y-  1) 
(y  -  xT 


(4.13) 
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Next  consider  a  purely  hypothetical  reversible  device  which 
is  not  an  ejector  at  all  but  a  type  of  ideal  turbomachine.  This 
device  receives  two  streams  of  fluid  identical  in  all  respects  to 


the  streams  at  stations  1  and  2  of  the  ideal  ejector.  It  dis¬ 
charges  a  stream  identical  in  all  respects  to  the  stream  at  sta¬ 
tion  3  of  the  ideal  ejector  except  for  the  pressure  pR3  which 

★ 

is  higher  than  p3  •  Again  the  density  p  is  taken  as  constant. 
The  reversible  device  operates  adiabatically  and  hence  isentropi- 
cally.  Under  these  conditions  the  following  energy  equation  holds. 
Recall  that  P2  -  p^  • 


xm3  ^  +  jPV]_2)  +  (1  -  x)  m3  ^  +  |pv22) 

=  m3  (p3R  +  ^Pv32)  (4.14) 

The  above  equation  can  be  simplified  by  subtracting  from  it, 
term  by  term,  the  following  identity. 

xm3p1  +  (1  -  x)m3P1  *  m3P1  (4.15) 

Upon  carrying  out  the  above  subtraction,  then  dividing  through 

.  2 

by  the  reference  power  m3  (V3  /2)  and  simplifying,  we  readily 
find  the  pressure  rise  coefficient  for  the  reversible  device  to 
be 


lp3R  -  p3> 

rpV32 


'PR 


(1  -  x)V 

(y  -  x) 


1 


(4.16) 
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Comparison  of  this  result  with  Eq.  (4.13)  shows  that  even  the 

* 

ideal  ejector  has  a  certain  inherent  mixing  loss,  call  it  cm  , 
as  compared  with  a  truly  reversible  device.  Thus 


°m  “  (cPr  -  V 


(4.17) 


Of  course  the  mixing  loss  cm  of  the  actual  ejector  will 
b.  somewhat  larger  than  cm*  and  u  i3  advantageous  to  empress 
it  in  the  form 


'm  nT 


(4.18) 


where  nE  /  which  we  term  ejector  effectiveness,  is  a  dimension¬ 
less  parameter  smaller  than  unity  which  must  be  determined  from 
experimental  data. 

Eq.  (4.18)  identifies  one  of  the  component  losses  that  reduce 
the  actual  pressure  rise  coefficient  cp  below  its  ideal  rever¬ 
sible  value  cpR  .  Another  loss  is  that  associated  with  separa¬ 
tion,  friction  and  secondary  flow  effects  within  the  ejector;  let 
us  denote  it  by  symbol  cfE  .  We  have  seen  that  a  third  component 
is  simply  the  gross  electrical  power  output  as  expressed  by  the  co¬ 
efficient  cQ  .  Thus  the  overall  energy  balance  within  the  ejec¬ 
tor  requires  that 


cp  =  CPR  "  n_  (cPR  cp  }  cfE  "  ce 
£ 


(4.19) 


The  resulting  pressure  rise  through  the  ejector  as  expressed 
by  Eq.  (4.19)  must  be  just  great  enough  to  offset  the  losses  in 
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the  return  passage.  These  again  are  losses  associated  with  separa¬ 
tion,  friction  and  secondary  flows;  let  us  denote  them  by  coeffi¬ 
cient  CjR  .  Accordingly,  for  the  return  passage  we  may  write 

cp  =  cfR  (4.20) 

We  now  eliminate  Cp  between  Eqs.  (4.19)  and  (4.20)  and  bring 
all  terms  to  one  side.  We  also  introduce  the  following  definition 
for  the  total  friction  loss  in  the  fluid  flywheel,  namely, 

cf  =  cfE  +  cfR  (4*21) 

The  result  is 

CPR  “  nE  (<=PR  cp  )  +  cf  “  Ce  =  0  (4.22) 

Now  using  Eqs.  (4.13),  (4.16)  and  (3.12)  to  express  the  above 
components  in  terms  of  variables  x  and  y  ,  we  obtain 


This  is  the  overall  energy  balance  equation  for  the  "fluid 
flywheel"  and  it  is  fundamental.  Notice  that  the  term  in  the  first 
pair  of  curly  brackets  represents  the  ideal  power  output  of  a  re¬ 
versible  device,  the  term  in  the  second  pair  of  curly  brackets 
represents  the  sum  of  the  fluid  power  losses  caused  by  mixing  and 
friction  and  the  term  in  the  third  pair  of  curly  brackets  represents 
the  gross  electrical  power  output. 
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Eq.  (4.23)  suggests  the  utility  of  defining  a  "conversion 


efficiency"  ncv  as  the  ratio  of  the  gross  electrical  power  cg 
to  the  ideal  power  cpR  .  Thus 


V  xy2 

(i  +  2xr 

xy2  ♦  &  - 

(y  -  xr 


(4.24) 


Notice  that  low  values  of  the  parameter  v  ,  which  is  ultimately 
determined  by  the  electrical  breakdown  limit,  will  necessarily 
produce  low  conversion  efficiencies. 

In  connection  with  Eq.  (4.23)  and  (4.24),  notice  that  x 
always  falls  between  zero  and  unity  while  y  is  always  greater 
than  unity  but  has  no  preassigned  upper  limit. 

Assuming  that  parameters  nE  ,  and  v  are  known  or  can 

be  estimated,  Eq.  (4.23)  then  fixes  a  corresponding  unique  relation 
between  variables  x  and  y  .  This  may  be  described  by  a  line  in 
the  xy  plane.  All  subsequent  calculations  must  be  confined  to 
points  that  lie  along  this  line. 

A 

If  the  area  ratio  __1  be  already  fixed,  then  we  have  from 

A3 

Eq.  (4.8)  the  additional  constraint 


Of  course  for  a  fixed  value  of  (Aj/A^) ,  Eq.  (4.25)  also 
represents  a  line  in  the  xy  plane.  In  fact  it  is  a  straight 


line  of  slope  (A^/A^)  which  passes  through  the  origin.  The  one 
or  more  points  of  intersection  of  this  straight  line  with  the 
curve  associated  with  Eq  (4.23)  fix  the  conditions  at  which  steady 
operation  of  the  fluid  flywheel  is  possible.  The  location  of  such 
a  point  determines  the  values  of  x  and  y  and  hence  ultimately 
determines  all  other  performance  characteristics  of  the  system  as 
well. 

Notice  that  Eq.  (4.23)  lends  itself  to  further  algebraic  sim¬ 
plification  which  makes  the  subsequent  numerical  solution  easier. 
On  the  other  hand  the  terms  in  the  simplified  result  cannot  be 
given  the  clear  physical  interpretation  that  is  possible  with 
Eq.  (4.23)  itself.  After  simplification,  Eq.  (4.23)  becomes 


(1  +  2X)‘ 

For  purposes  of  numerical  solution,  it  is  useful  to  cast 
Eq.  (4.26)  into  a  more  tractable  form  as  follows.  Let 


(4.27) 


(4.28) 


(4.29) 
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(4.30) 


With  this  notation,  Eq.  (4.26)  may  be  rewritten  as  follows, 
where  F  is  a  function  that  must  equal  zero.  Thus 


uaii' 

(y-x)2 


(l+5x) 


(l+2x) 


:Tf  *  A4  " 


(4.31) 


2  2 

Next  we  multiply  through  by  (1  +  2x)  (y  -  x)  and  denote 
the  resulting  function  by  G.  Upon  expanding  the  right  side, 
regrouping  terms  and  simplifying,  we  find  that  G  is  a  quartic 
in  x.  The  coefficients  of  this  polynomial  turn  out  to  be 


B4  =  -  4A1-4(3-2y)A2-5A3-4A4  (4.32) 
B3  =  4 (y-l)A1-4y(y-2)A2+(10y-l)A3+4 (2y-l)A4  (4.33) 
B2  =  (4y-l)A1-(4y2-2y-5)A2-y(5y-2)A3-(4y2-8y+l)A4  (4.34) 
B1  *  yA1-(y2+DA2-y2A3-2y(2y-l)A4  (4.35) 

B0  =  "A2‘y2A4  (4*36) 
The  function  G  may  then  be  written 


G  =  B4x4+  B3x3+B2x2+B1x+B0  =  0  (4.37) 

The  unknown  root  x  can  be  found  by  Newton's  method.  For 
this  purpose  we  require  also  the  derivative  of  G  which  is 


=  4B 


,x3+3B,x2+2B-x+B1 


(4.38) 
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Let  x 


be  any  trial  value  of  x  and  let  and  G' 

n  ■*  n  n 

be  the  corresponding  values  computed  from  Eqs.  (4.37)  and  (4.38). 
Then  according  to  Newton’s  method,  the  next  trial  value  xn+^ 
should  be  taken  as 


(4.39) 


This  procedure  converges  rapidly  to  the  true  root  x  . 

This  method  was  used  in  calculating  the  numerical  example  which 
is  presented  in  a  later  section. 
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>.  Thermodynamic  Performance 

As  shown  in  the  schematic  temperature  entropy  diagram.  Fig 
5 ;  1 ,  the  tluid  flywheel  receives  a  supply  m^  (kg/sec)  of  super¬ 
heated  steam  at  condition  p  ,  T  and  discharges  saturated  or 

o  o 

nearly  saturated  liquid  at  condition  p4  ,  T4  .  It  gives  up  heat 
to  the  ambient  atmosphere  at  absolute  temperature  TA  and  delivers 
a  gross  electrical  power  output  P&  (watts) .  It  also  requires 
an  input  of  electric  power,  call  it  ,  to  energize  the  injec¬ 

tor.  This  is  normally  a  small  percentage  of  the  gross  electrical 
output  and  is  conveniently  expressed  by  neans  of  an  excitation 
efficiency  nx  such  that 


(watts) 


(5.1) 


(continued  on  page  5-3) 


After  being  condensed  and  leaving  the  fluid  flywheel,  the 
liquid  water  must  be  pumped  back  into  the  boiler.  For  an  ideal 
pump  of  100%  efficiency  this  would  require  the  expenditure  of  shaft 
power  in  the  amount 

Pp*  *  m^  (pQ  -  p4)  ^  (watts)  (5.2) 


where  v4  denotes  the  specific  volume  of  the  liquid  at  state  4. 
For  a  real  pump  of  efficiency  n_  the  corresponding  power  is 


P 


P 


(watts) 


(5.3) 


Thus  it  is  seen  that  for  the  overall  cycle, 
power  output  may  be  written  as 

pe"  *  (nxPe  -  pp*/op)  (watts) 


the  net  useful 


(5.4) 


The  rate  of  heat  input  to  the  steam  passing  through  the  boiler 
and  superheater  is 

Q  3  n»i  (hQ-h4)  -  Pp  /np  (watts)  (5.5) 

where  hQ  and  h4  denote  the  enthalpies  (joule/kg)  at  states 
o  and  4 . 

For  the  hypothetical  case  of  100%  pump  efficiency,  we  denote  the 

•  ★ 

corresponding  heat  rate  calculated  from  Eq.  (5.5)  by  symbol  Q  . 

This  differs  only  very  slightly  from  Q  . 
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The  overall  cycle  efficiency  now  becomes 


(5.6) 


It  is  of  interest  to  compare  the  above  performance  of  the 
actual  EGD  device  with  the  ideal  performance  of  a  hypothetical 
reversible  device  operating  under  the  same  input,  output  and  am¬ 
bient  conditions.  According  to  the  theory  of  availability  in 
steady  flow,  the  maximum  electrical  power  that  could  be  produced 
by  such  an  ideal  device  is 

pe  *  ^  | (hQ  -  h4)  -  TA  (sQ  -  s4)j  (watts)  (5.7) 

where  sQ  and  '  s4  denote  the  entropies  at  states  o  and  4, 
respectively  (joule/kg  °K) . 

The  cycle  efficiency  of  this  ideal  device  would  be 

★ 

it  1  *  .  p  \ 

nc  =  (Pe  p  (5*8) 

C  Q  e 

It  is  also  instructive  to  compare  the  net  output  of  the  real 
system  with  that  of  the  above  ideal  and  reversible  device.  We 
term  this  ratio  the  relative  output.  Thus 


nrel 


(5.9) 


All  of  the  foregoing  relations  can  be  put  into  dimensionless 
form  by  dividing  all  quantities  having  the  units  of  power  by  a 
suitable  reference  power.  In  the  previous  analysis  of  the  fluid 


flywheel,  it  was  natural  and  convenient  to  choose  this  reference 
as  the  kinetic  power  of  the  total  stream  at  the  ejector  exit, 
station  3.  While  it  would  be  possible  to  retain  this  same  reference 
in  this  section,  it  is  advantageous  to  choose  here  as  reference, 

the  kinetic  power  of  the  primary  jet  at  the  ejector  inlet,  station 

% 

1.  This  choice  reduces  the  final  results  to  a  simpler  form  and 
makes  them  somewhat  easier  to  interpret.  Thus  the  reference  power 
is  now 


In  computing  this  reference  power  it  is  permissible  to  assume 
that  state  1  is  fixed  by  the  two  conditions 


S1  *  so 


(5.n; 


Pi  -  p4 


(5.12) 


In  connection  with  Eq.  (5.10)  it  is  also  instructive  to  note 


that 


(hQ  -  hx)  = 


Vl2  I  *o\2 


i  ♦  w 


(5.13) 


Recall  also  that  if  we  introduce  the  optimizing  condition 
developed  earlier,  namely, 


2  2  1 
M  =  M  .  =i 
1  lopt  2 


(5.14) 
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then  Eq  (5.13)  simplifies  further  to 


(h. 


-V-11 


~o 

Ty+3T 


(5.15) 


Incidentally,  it  is  also  worth  noting  that  the  overall  pres¬ 
sure  ratio  across  the  fluid  flywheel  is 


Moreover,  if  we  observe  the  restriction  of  Eq  (5.16)  this  becomes 


(5.17) 


For  example,  for  y  *  1.3  this  gives 

©O"  °'731 

This  result  expresses  one  of  the  fundamental  limitations  of 
the  fluid  flywheel  concept,  namely  the  restriction  to  a  fixed 
pressure  ratio.  The  reason  for  this  restriction  was  explained 
earlier  in  connection  with  Fig.  2.1. 
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Upon  simplifying  and  reducing  the  various  foregoing  thermo¬ 
dynamic  relations  in  the  manner  explained  above,  we  obtain  the 


following  results. 


riii  (hQ-h1) 


,  .  (W  (so-S4> 
1  +  — — 


Reversible 
Electric  Power 


m.  (h  -h. ) 
1  o  1 


=  K 


•  * 
Q 


m.  (h  -h. ) 
1  o  1 


=  K„ 


m.  (h  -h, ) 
1  o  1 


=  K, 


(Pp~P4)v4 
(h o~hl] 


(W 


-  K. 


<VV 


nP 


Reversible 
Pump  Power 

Reversible 
Heat  Input 

Actual  Heat 
Input 


e 


=  K 


V 


(l+5x) 
(1+2  x)  2 


Gross  Electric 
Power 


ll 


(h  -h. ) 
1  o  1 


Net  Electric 
Power 

Cycle  Efficiency 


i  *  * 

7*  (Ke  -  Kp  > 


rel 


n  K  -  K  /n  , 
x  e  p  p ' 


K  _  K 
e  Ap 


Reversible 
Cycle  Efficiency 

Relative  Output 


5.19) 

5.20) 

(5.21) 

5.22) 

5.23) 

5.24) 

5.25) 

5.26) 

5.27) 


5-7 


Once  Eq.  (4.37)  of  the  previous  section  is  solved,  it  becomes 
a  relatively  straight  forward  matter  to  apply  the  relations  of 
the  present  section  to  calculate  the  various  thermodynamic  per¬ 
formance  parameters.  Recall  that  Eq.  (4.37)  defines  a  whole  series 
of  points  in  the  xy  plane.  By  calculating  the  thermodynamic 
performance  parameters  of  this  section  for  such  a  series  of  points, 
and  by  comparing  the  respective  results,  it  becomes  possible  to 
choose  that  optimum  design  point  which  yields  the  best  overall 
trade  off  among  the  various  desired  performance  characteristics. 

So  far  only  an  initial  trial  calculation  of  this  type  has 
been  made  on  the  basis  of  a  plausible  initial  set  of  key  para¬ 
meters.  This  result  is  summarized  in  a  later  section. 
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6 .  Summary  of  Principal  Calculation  Formulas 


The  principal  equations  developed  in  the  preceding  sections 
are  summarized  below  for  easy  reference.  They  are  listed  without 
further  explanation  in  the  approximate  order  in  which  they  would 
be  used  in  a  typical  calculation. 


Input  Parameters 

1 


(hQ-h1)  (where  s1  =  sQ,  px  =  p4) 


(W 


(P0-P4)v4 


<W 

TVE7T-  kp 

o  i 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 

(6.7) 

(6.8) 
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(6.9) 


Solution  of  Energy  Balance  Equation 


A1 

.  2(y-i>2 

(6.11) 

A2 

■( 

•H’’ 

(6.12) 

A3 

3 

vy2 

(6.13) 

A4 

- 

(°<  ■  ^  *  ‘) 

(6.14) 

B4 

= 

-4A^-4  (3-2y)  A2-5A3~4A4 

(6.15) 

B3 

s 

4 (y-l)A1-4y (y-2)A2+(10y-l) A3+4 (2y-l)A4 

(6.16) 

B2 

= 

(4y-l)A1- (4y2-2y-5)A2-y (5y-2) A3- (4y2-8y+l)A4 

(6.17) 

B1 

s 

yA1- (y2+l) A2-y2A3-2y (2y-l) A4 

(6.18) 

Bo 

= 

-A2-y2A4 

(6.19) 

G 

- 

4+Bjx  3+B2X  2+Bxx+Bo 

(6.20) 

G' 

s 

4B4x3+3B3x2+2B2x+B1 

(6.21) 

xn+l 

3 

Gn 

xn  -  £—i  (Iterate  until  x  converges) 

(6.22) 

I 

L 
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Output  Parameters 


* 


II 

Ke 

(l+5x) 

i  nxv  2 

1  (l+2x) i 

-!ej 

nP 

|  Net  Electric 

|  Power 

(6.23) 

nc 

s 

K 

e 

kq 

Cycle  Efficiency 

(6.24) 

in 

K  " 
e 

Relative  Output 

(6.25) 

nrel 

★  W 

<Ke  ‘Kp  > 

Ke 

s 

1  (1+5X) 

1  (l+2x)2j 

Gross  Electric 

Power 

(6.26) 

kpr 

! 

\l  +  iiisil. 

(  x  f y-x ) z 

xy2  1 

Ideal  Pressure 

J  Rise 

(6.27) 

cv 

= 

^e 

kpr 

Conversion 

Efficiency 

(6.28) 

Note:  Power  coefficients  denoted  by  symbol  K  are  expressed 

in  terms  of  the  kinetic  power  at  station  i.  Those  denoted  by 
symbol  c  are  expressed  in  terms  of  the  kinetic  power  at  station 
3.  The  conversion  between  these  two  forms  is  simply 

K  =  -—s’  ( 6 . 29 ) 

xy 
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Then  from  Eq.  (6.3)  we  obtain 


In  the  absence  of  adequate  test  date,  the  best  estimate  of 
u  that  is  possible  at  this  time  is  simply 

U  «  1  (7.4) 

According  to  Eq.  (6.4)  this  gives  the  estimate 

v  -  m|»6  -  (1)  (1.32)  (0.0238)  =  0.032  (7.5) 

The  value  of  the  flow  loss  factor  cf  for  the  annulus  may 
be  estimated  from  the  data  in  Ref.  (2)  on  losses  in  pipe  bends. 

At  a  Reynolds  number  of  200,000  or  above,  these  losses  are  mainly 
a  function  of  R/r,  where  R  is  bend  radius  and  r  is  pipe 
radius.  For  the  torus  in  Marks'  design  R/r  =  5.7.  Ref.  (2) 
shows  that  at  this  value  of  R/r  the  loss  factor  for  a  180° 
pipe  bend  is  about  0.4.  Thus  for  a  360°  turn  the  loss  factor 
would  presumably  be  twice  the  value  or  approximately 

c f  *  0.8  (7.6) 

There  are  qualitative  grounds  for  hoping  that  the  loss  in 
a  true  torus  might  be  somewhat  smaller  than  the  above  pipe  bend 
data  would  suggest,  but  there  is  no  hard  data  available  on  this 
point.  Thus  Eq.  (7.6)  represents  a  plausible  and  conservative 


working  hypothesis  that  can  be  advanced  at  this  time. 

Values  of  nE  can  be  estimated  from  available  test  data 
on  ejector  performance.  See,  for  example.  Ref.  (3) .  These 
data  have  net  yet  been  reviewed  critically.  All  we  can  say  at 
present  on  the  basis  of  general  engineering  judgment  and  experi¬ 
ence  is  that  for  the  range  of  ejector  design  parameters  of  in¬ 
terest  in  the  present  context,  nE  is  unlikely  to  lie  above 
0.9  or  below  0.5.  It  is  thought  reasonable  as  a  tentative  ini¬ 
tial  working  hypothesis  to  choose 

nE  «  0.8  (7.7) 

Also  required  for  the  calculation  are  certain  thermodynamic 
properties.  These  are  summarized  in  Table  7.1.  These  data  hap¬ 
pen  to  be  in  English  units  but  this  is  permissible  because  the 
final  results,  being  dimensionless,  are  independent  of  the  system 
of  units  used.  For  some  purposes  the  steam  was  treated  approxi¬ 
mately  as  a  perfect  gas  with  a  ratio  of  specific  heats  y  =  1.3. 

The  information  in  Table  7.1,  along  with  the  formulas  of 
Section  6,  permit  the  various  key  input  parameters  to  be  esti¬ 
mated  as  summarized  in  Table  7.2.  These  values  are  regarded  as 
the  best  estimates  that  are  possible  at  the  present  time.  These 
input  data  suffice  to  permit  a  corresponding  performance  calcula¬ 
tion  to  be  carried  out.  This  has  been  done  and  the  main  results 
are  summarized  in  Table  7.3  and  in  Fig.  7.1.  Notice  that  each 
row  or  Table  7.3  corresponds  to  a  different  value  of  the  ejector 
area  ratio  A^/A^. 


Table  7.1  Thermodynamic  Properties 


State 

0 

1 

4 

Ambient 

P 

psia 

1500 

1100 

1100 

— 

T 

°R 

1063.9 

1016 

1016 

530 

h 

Btu/lbm 

1180.7 

1156.3 

557.4 

— 

s 

Btu/lbm°R 

1.3470 

1.3470 

0.7575 

— 

P 

lbm/ft3 

— 

— 

45.45 

Table  7.2  Estimated  Values  of  Key  Input  Parameters 


Y 

= 

1.3 

np 

= 

0.9 

U 

= 

1 

* 

K 

e 

= 

12.7 

V 

= 

0.032 

* 

KP 

= 

0.067 

nE 

= 

0.8 

nr 

Ko 

= 

25.5 

Cf 

= 

0.8 

kq 

= 

25.5 

s 

1 

Hr 

— 

0.50 

Table  7.3  Estimated  Performance  of  EHD  Generator 
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The  curve  ABCDEF  shown  in  Fig.  7.1  is  the  locus  of  points 
x,y  which  satisfy  the  basic  energy  balance  relation  of  the 
ejector  and  fluid  flywheel/  that  is,  Eq.  (4.23).  In  view  of 
the  tentative  nature  of  the  values  assigned  to  parameters  v, 
nE  and  Cf,  the  accuracy  of  this  curve  is  somewhat  uncertain, 
but  its  general  qualitative  features  are  undoubtedly  correct. 

Any  straight  line  througn  the  origin,  such  as  line  OEB,  for 
example,  represents  Eq.  (4.25),  and  the  inverse  of  its  slope 
represents  the  area  ratio  A^/A 3  .  Notice  that  such  a  straight 
line  may  in  general  intersect  the  curve  at  two  points  such  as 
E  and  B  ,  for  example.  It  can  be  shown  that  only  the  lower 
point  E  represents  stable  equilibrium;  point  B  represents 
unstable  equilibrium  and  should  be  disregarded.  For  this  reason 
the  entire  dashed  portion  ABC  of  the  curve  must  also  be  dis¬ 
carded.  Only  the  solid  portion  CDEF  has  any  practical  signifi¬ 
cance  . 

Notice  that  the  tangent  point  C  represents  the  lowest  value 

of  aj_/A3  and  the  highest  value  of  y  that  can  be  achieved  in 

steady  state  operation.  Recall  that  the  gross  electric  power  is 

2 

proportional  to  the  factor  (l+5x) / (l+2x)  .  It  can  be  shown  that 
this  factor  reaches  its  peak  value  of  1.04  at  x  *  0.1  .  It  is 
therefore  desirable  to  choose  the  operating  point  on  the  curve 
as  close  to  x  =  0.1  as  possible.  In  this  case  the  nearest  point 
is  at  D  where  x  =  0.36.  Note  that  this  is  the  lowest  value  of 
x  that  can  be  attained  for  the  assumed  values  of  v ,  nE  and  cf  . 
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The  corresponding  value  of  the  factor  (l+5x)  /  (l+2x)  ^  is  0.95. 

Notice  from  the  data  in  Table  7.2  that  an  ideal  reversible 
device  operating  on  the  given  basic  thermodynamic  cycle  would 
produce  a  dimensionless  net  power  output  of  amount 

(Kg*  -  Kp*)  =  12.7  -  0.07*  12.6  (7>g) 

> 

and  would  produce  an  ideal  cycle  efficiency  of  amount 

nc*  =  0.50  (7.9) 

These  theoretical  performance  limits  of  the  ideal  reversible 
device  are  in  themselves  very  favorable.  Unfortunately,  the  per¬ 
formance  of  the  actual  EHD  system  as  summarized  in  Table  7.3  is 
seen  to  be  very  unsatisfactory.  In  fact  the  net  power  outputs 
and  cycle  efficiencies  are  negative  for  all  possible  operating 
points!  This  means  that  even  at  100  atmospheres  pressure,  the 
phenomenon  of  electrical  breakdown  so  severely  limits  the  gross 
electrical  power  output  that  it  it  insufficient  even  to  meet  the 
demand  for  power  to  drive  the  pump,  let  alone  provide  any  useful 
output. 

It  might  be  supposed  that  a  radical  increase  in  the  electrical 
breakdown  strength,  if  it  could  be  achieved,  would  in  itself  trans¬ 
form  this  pessimistic  picture  into  an  optimistic  one,  but  this  is 
not  entirely  the  case,  at  least  not  for  the  present  EHD  scheme. 
Suppose  for  example,  that  the  breakdown  parameter  v  could  be  in¬ 
creased  by  a  factor  of  roughly  30  so  as  to  bring  the  gross  elec¬ 
trical  power  coefficient  up  say  to 

Ke  *»  1  (7.  10) 
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In  that  case  the  cycle  efficiency  would  still  be  only  about 


n 


c 


(1  -  0.07)  *  0.037 


(7.11) 


or  less  than  4%.  The  main  reason  for  this  unsatisfactory  result 
is  that  the  present  cycle  squanders  available  energy  by  dumping 
heat  from  the  condenser  at  high  temperature  to  the  ambient  air 
at  low  temperature  without  extracting  any  useful  power  from  it. 

The  irreversibilities  inherent  in  the  ejector  also  waste  available 
energy.  Moreover,  friction  losses  in  the  fluid  flywheel  appear 
to  be  disappointingly  high.  As  a  reflection  of  these  losses, 
note  from  Table  7.3  that  the  conversion  efficiency  ncv  remains 
below  4%  for  all  possible  stable  operating  points. 

It  is  also  worthwhile  at  this  point  to  review  the  various 
factors  that  affect  the  important  parameter  v  which  so  severely 
limits  the  electrical  output.  These  factors  are  revealed  clearly 
in  Eqs .  (6.2),  (6.3)  and  (6.4).  It  is  convenient  to  combine  them 
here  in  the  following  way,  assuming  a  -  0.  Thus 


!Y+1 
v  [yTj] 

As  y  varies  from  1  to  1.67,  the  function  of  y  in  the 
curly  brackets  above  varies  from  1.82  to  0.97.  Hence  there  is 
no  possibility  of  radically  changing  the  order  of  magnitude  of 
v  by  manipulating  y  . 


f£CB2p0 


(7. 12) 
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Similar  considerations  apply  to  the  quantity  e  .  It  has 
been  suggested,  see  Re.  4,  that  e  might  be  increased  by  using 
a  polarizable  aerosol.  But  the  overall  increase  that  can  be 
achieved  in  this  way  is  expected  to  be  less  than  a  factor  of  ten 
and,  as  we  have  seen,  this  in  itself  is  not  enough  to  overcome 
fully  the  limitations  of  the  present  EHD  scheme. 

An  obvious  way  to  increase  v  is  to  increase  pQ  .  We  have 
noted,  however,  that  even  at  100  atmospheres  pressure,  v  is 
still  much  too  small.  Any  increase  beyond  this  value  is  not 
feasible  for  two  reasons.  Firstly,  structural  and  practical 
problems  become  troublesome  at  such  high  pressures.  Secondly 
the  linear  relation  between  electrical  breakdown  strength  and 
fluid  density  as  given  by  Eq.  (2.5)  breaks  down  near  pressures 
of  100  atmospheres.  See  Ref.  1. 

The  two  remaining  factors  at  our  disposal  are  CB  and  TQ  . 

2 

Note  that  v  is  proportional  to  the  quadractic  factor  (CB/TQ)  . 

The  breakdown  constant  C_  is  a  property  of  the  medium.  The 
question  whether  any  aerosol  mixture  exists  that  has  a  substan¬ 
tially  higher  value  of  CB  has  not  been  investigated.  It  is 
recommended  that  such  an  investigation  be  carried  out. 

Finally  consider  the  factor  TQ  .  In  the  present  scheme, 
charged  particles  of  low  mobility  are  created  by  condensation. 
This  requires  that  TQ  be  high  enough  to  produce  about  5% 
liquid  at  state  1.  This  restriction  on  TQ  could  be  eliminated 
by  abandoning  the  condensation  method  and  utilizing  instead  the 
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direct  injection  of  liquid  droplets  or  solid  particles.  This 
change  also  entails  changing  the  carrier  fluid  from  steam  to 
air  or  some  other  noncondensing  gas .  Note  that  a  decrease  in 
T  Q  from  1016  °R  to  530  °R  would  increase  v  by  a  factor  of 
about  3.7. 

It  seems  possible  that  some  combination  of  the  above  methods 
could  perhaps  increase  v  by  a  factor  of  say  20  or  so.  We  have 
seen  that  this  in  itself  would  not  make  the  present  fluid  fly¬ 
wheel  configuration  successful.  Nevertheless,  an  improvement  of 
this  order  might  very  well  suffice  to  make  some  alternative  EHD 
scheme  practical. 
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8.  Conclusions  and  Recommendations. 


Broadly  speaking,  one  of  the  inherent  limitations  of  EHD 
power  generation  is  the  severe  restriction  imposed  by  the 
phenomenon  of  electrical  breakdown.  Because  of  this  effect, 
the  electrical  work  that  can  be  extracted  per  unit  mass  of  work¬ 
ing  fluid  is  very  small.  In  general  there  are  two  distinct  ap¬ 
proaches  for  attempting  to  overcome  this  limitation.  One  method 
is  based  on  multi-staging,  that  is,  on  passing  the  same  working 
fluid  through  a  whole  series  of  successive  electrical  conversion 
sections.  The  other  method  is  to  adopt  an  ejector  system  which 
utilizes  a  primary  stream  of  high  velocity  but  low  mass  flow 
rate  to  drive  a  secondary  stream  of  low  velocity  but  high  mass 
flow  rate.  The  conversion  of  flow  energy  to  electric  work  is 
accomplished  in  the  resulting  low  velocity  stream.  Although 
electrical  work  per  unit  mass  of  low  velocity  fluid  remains  low, 
the  electrical  work  per  unit  mass  of  high  velocity  primary 
fluid  may  be  boosted  to  a  sufficiently  high  value.  The  effective¬ 
ness  of  this  scheme  is  limited  to  some  extent  by  the  character¬ 
istic  ejector  mixing  losses. 

In  this  report,  a  proposed  type  of  single  stage  EHD  genera¬ 
tor  is  analyzed  which  utilized  the  ejector  principle.  However, 
this  design  combines  the  ejector  and  the  electrical  conversion 
section  into  a  single  unit;  in  particular,  charged  particles  are 
introduced  into  the  primary  stream  only  and  are  subsequently 
collected  at  the  ejector  exit  section. 

The  analysis  shows  that,  unfortunately,  this  type  of  design 
gives  extremely  poor  performance.  In  fact  the  electrical  power 
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generated  is  insufficient  to  drive  the  feedwater  pump,  let  alone 
produce  any  net  useful  output.  This  poor  result  can  be  attribu¬ 
ted  to  the  fact  that  the  charged  particles  are  introduced  into 
the  primary  stream,  thus  generating  a  very  intense  local 
electrical  field  at  this  point. 

An  obvious  remedy  for  this  difficulty  is  to  place  the  elec¬ 
trical  conversion  section  after  the  ejector  where  the  electrical 
field  can  be  much  lower.  The  Appendix  evaluates  this  revised 
design.  It  is  shown  in  the  Appendix  that  under  the  most  opti¬ 
mistic  assumptions,  the  cycle  efficiency  is  less  than  4%  which 
may  be  insufficient  to  compete  with  alternative  power  systems. 

It  may  be  possible  to  multi-stage  the  fluid  flywheel 
generators  and  thereby  show  an  order  of  magnitude  improvement 
in  performance.  It  is  recommended  that  systematic  parametric 
studies  be  carried  out  to  determine  the  performance  possi¬ 
bilities  of  such  a  multi-stage  system  and  to  estimate  optimum 
design  values  of  the  key  parameters.  It  would  also  be  ad¬ 
visable  in  the  interest  of  orderly  development  to  restrict  this 
proposed  next  step  to  the  same  basic  medium,  namely  steam,  and 
to  the  same  basic  Rankine  cycle  as  used  in  the  previous  studies. 
Of  course  there  is  no  assurance  that  a  multi-stage  design  will 
necessarily  prove  successful  because,  even  if  the  theoretical 
cycle  performance  is  favorable,  practical  considerations  of 
bulk  and  complexity  might  nullify  the  theoretical  advantages. 
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APPENDIX:  Performance  Characteristics  of  Revised  EHD  Generator 


The  analysis  in  the  main  text  clearly  indicates  that  intro¬ 
ducing  charges  in  the  primary  stream  only,  limits  the  thermody¬ 
namic  state  1  to  a  relatively  high  fixed  pressure  ratio 
(p^/pQ  =  0.731  for  =  0.71).  An  obvious  improvement  is  to 
uncouple  the  ejector  from  the  electrical  conversion  section. 

This  is  done  in  the  following  analysis  in  which  the  electrical 
conversion  section  is  located  after  the  ejector. 

The  work  in  this  Appendix  is  based  on  a  fixed  upper  pressure 
pQ  =  1500  psia  and  a  fixed  steam  quality  of  0.950  at  station  1 
(See  Fig.  5.1).  However,  the  back  pressure  p^  is  allowed  to  take 
on  a  large  number  of  possible  values  ranging  from  1400  psia 
down  to  14.7  psia.  At  any  given  back  pressure,  the  area  ratio 
A^/A^  is  allowed  to  vary  as  required  to  obtain  best  performance. 
It  is  found  that  best  performance  corresponds  to  the  minimum 
possible  value  of  the  mass  flow  ratio  (ift^/ih^)  =  x.  This  is 
governed  by  the  energy  balance  equation  for  the  fluid  flywheel. 

Three  cases  are  considered,  as  listed  in  more  detail  on 
page  2A.  Case  A  represents  the  present  best  estimate  of  the 
three  key  parameters  c^  and  CQ.  Case  B  represents  a  much 

more  optimistic  estimate  of  parameters  nR  and  c^,  perhaps  values 
which  can  never  be  achieved  practically.  Case  C  retains  these 
optimistic  values  of  nR  and  Cf  and  assumes  in  addition  that  the 
electrical  breakdown  constant  CR  can  be  increased  by  a  factor 
of  /To.  Thus,  case  C  represents  optimism  plus  a  fundamental 
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REVISED  EHD  ANALYSIS 


BASIS 

*  Design  change:  electrical  power  section  now  placed  after 
ejector. 

*  Rankine  cycle  with  - 

pQ  -  1500  psia 

steam  quality  at  state  1  -  0.95 


OPTIMIZING  PARAMETERS 


*  Pressure  p^  in  fluid  flywheel 

*  Ejector  area  ratio  A^/A^ 

CASES  STUDIED 

ABC 


Ejector  Effectiveness  nE 

Friction  Factor  cf 

Electrical  Breakdown  Cg 

,  Constant 
(nr  °K/Coulomb) 


0.8  0.9  0.9 

0.8  0.1  0.1 

9 . 49xl03  9 . 49xl03  /l0x9 . 49xl03 


scientific  "break-through" . 


The  single  most  important  result  of  the  analysis  is  the 
curve  of  cycle  efficiency  n  versus  back  pressure  p,  ,  Fig.  1. 
This  result  confirms  that  the  proposed  design  revision  has  indeed 
improved  the  performance  over  that  of  the  original  version.  Re¬ 
call  that  the  original  design  yielded  a  negative  cycle  efficiency! 
The  revised  design,  case  A,  product.-  positive  efficiencies  for 
values  of  p^  above  about  900  psia.  The  peak  efficiency  occurs 
approximately  at  p^  *  1400  psia.  Unfortunately,  this  peak  value 
is  still  much  too  low,  below  0.003.  As  expected,  case  B  shows  a 
marked  improvement  over  case  A,  but  nevertheless  nc  still  re¬ 
mains  below  0.01.  Even  the  extremely  optimistic  case  C  gives  a 
peak  cycle  efficiency  of  less  than  0.04. 

The  efficiency  curve  for  case  C  has  an  interesting  bimodal 
form  with  the  solution  above  p1  =  1000  psia  quite  different  from 
that  below  this  limit.  The  analysis  shows  that  the  ejector  cannot 
function  for  case  C  for  values  of  p^  above  about  1000  psia,  so 
the  solution  over  this  region  pertains  to  a  simple  EHD  conversion 
in  a  single  stream;  the  ejector  principle  cannot  be  utilized  in 
this  region.  It  is  ironic  that  the  highest  efficiency  shown, 
nearly  0.04,  occurs  at  about  p^  =  1400  psia  where  the  fluid  fly¬ 
wheel  concept  does  not  even  apply!  A  second  maximum  cycle  effi¬ 
ciency  of  about  0.026  occurs  at  p^  =  270  psia.  This  point  repre¬ 
sents  normal  operation  which  utilizes  the  fluid  flywheel  principle 
to  good  advantage . 


3A 


The  results  show  that  the  final  cycle  efficiency  that  can 


be  acheived,  and  the  optimum  back  pressure  p^  and  ejector  area 
ratio  (A^/A^)  that  will  produce  this  optimum  depend  critically 
upon  the  three  key  parameters  nE,  cf  and  Cg.  If  these  para¬ 
meters  can  be  estimated  accurately  then  design  optimization  and 
resulting  performance  can  be  accurately  predicted  and  controlled, 
otherwise  not. 

The  results  also  show,  however,  that  even  under  the  most 
optimistic  assumptions  regarding  nE,  c^  and  Cg,  the  cycle 
efficiency  is  disappointingly  low  -  less  than  0.041  It  is 
therefore  clear  that,  viewed  as  a  single  stage  device,  the  pro¬ 
posed  EHD  generator,  even  with  the  noted  design  improvements, 
has  performance  limitations  which  are  just  too  severe  to  permit 
it  to  compete  successfully  with  alternative  power  systems. 

On  the  other  hand,  the  possibility  still  remains  that  a 
multi-stage  version  of  this  generator  might  show  a  further  order 
of  magnitude  improvement  in  performance.  It  is  not  clear  at 
this  time  whether  the  ejector  and  fluid  flywheel  concepts  can 
be  incorporated  to  advantage  in  such  a  multi-stage  design. 

The  remainder  of  this  Appendix  summarizes  the  analysis  of 
the  revised  scheme.  The  notation  and  approach  is  generally 
similar  to  the  main  text.  The  derivations,  formulas,  tabula¬ 
tions  and  diagrams  are  complete  but  no  attempt  has  been  made 
to  integrate  them  into  a  smooth  flowing  text. 
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ELECTRIC  POWER  COEFFICIENTS 


P  2  V1 
5^  -  c  ■  vy  where  y=  («-) 

PK3  e  3 


=  K  =  —  where  x  =  (— ) 
PKI  6  X 


v  -  0.5837- 


!ho  -  V 


mi| 
m3 

(CASES  A  &  B) 


p  =  density,  lbm/ft 

(h  -  h. )  =  enthalpy  drop,  Btu/lbm 
o  1 


Note:  Above  constant  0.5837  is  based  on  a  value  of  the  break 

down  constant  of: 

C_,  =  9.49  x  10 3  m2°K/coulomb 

O 

An  assumed  increase  in  C„  by  a  factor  /T6’  would  increase  v  by 
a  factor  of  10.  This  would  give: 


v 


5.837 


P1 


hl} 


(CASE  C) 
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ENERGY  BALANCE 


OPERATING  LIMIT 

The  fluid  flywheel  can  function  only  if  parameters  nE#  cf 
and  v  satisfy  the  inequality 


cf  (-i  -  1)  +  v  (-i  +  1  +  c-)<l 

f  Hr,  f 


(1) 


OPERATING  RANGE 


Let 


C2  *  %  -  1  +  ^ 


C1  =  _1 
nE 


C  =  (-=■  +  1  +  cf) 

°  nE  f 


(2) 

(3) 

(4) 


R2  =  C?  -  4C-C 
1  ■do 


(5) 


A  solution  exists  for  any  valve  of  y  between  the  limits  - 


^min  ~  ^  2C 


C.  -  R 
J-  \ 


(6) 


=  ( 


C.  +  R 

*•  \ 


rmax  2C 


(7) 


MASS  RATIO  x: 


The  solution  for  mass  ratio  x  is  found  by  the  following 
sequence  of  calculations. 


=  2J3LZ..1). 


A-  = 


(“  "  1)y2 
nE 


A3  =  (-  -i-  +  1  +  cf  +  vy2) 


(8) 

(9) 

(10) 
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(11) 


B2  =  A1  +  (3  "  2y)A2  +  A3 

B1  =  yAl  "  (y2  “  3)A2  +  2yA3 
Bo  =  A2  +  y2A3 


(12) 

(13) 


S2  =*  B?  -  4B0B  (14) 

1  -4  0 

X  =  {\-i)  (151 

Best  performance  corresponds  to  the  minimum  possible  value 
This  is  best  found  by  numerical  trial  and  error. 


THERMODYNAMIC  PARAMETERS 


Let  p  =  power,  watts 

PK1  =  kinetic  power  at  station  1,  watts 


Then 


5^-  =  power  coefficient 
PK2 


*  -  ta)  (So  -  s4, 


Ke  -  a  + 


(ho  -  V 


„*  <Po  -  ?4)V4 

p  nr  -  h,  i 


o  "1 
(h  -  h,) 


KQ  =  {(h°  -  h  ) 
o  1 


A  J  *\ 


Reversible  Electric  Power 

Reversible  Pump  Power 

Heat  Input 


Pe 


m. 


n  v  k* 

{ir  '  (ho  "  V 


Specific  Output  Watt  Sec/Kg 


n  v  K  i 
n  =  {— - E}_I 

C  x  nP  KQ 


Cycle  Efficiency 


PRINCIPAL  PERFORMANCE  PARAMETERS 
Values  used  in  calculations- 

nxv 

Excitation  Factor  —r—  =  1.00 


Pump  Efficiency 


n  =  0.90 

r' 
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THERMODYNAMIC  PROPERTIES 
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THERMODYNAMIC  COEFFICIENTS 


Pi 

psia 

V 

K 

1400 

0.4063 

54.94 

1300 

0.1721 

26.53 

1200 

0.1001 

17.65 

1100 

0.06274 

12.75 

1000 

0.04244 

10.04 

900 

0.02958 

8.151 

800 

0.02073 

6.769 

700 

0.01455 

5.716 

600 

0.01010 

4.876 

500 

0.006839 

4.178 

400 

0.004430 

3.574 

300 

0.002663 

3.037 

200 

0.001391 

2.537 

100 

0.0005091 

2.016 

14.696 

0.00004396 

1.374 

Notes : 

Ah 

=  (hQ  -  hx) 

AH 

-  (h°  -  Vp, 

=  14.676 

kq 

KS 

(— ) 
'ah' 

110.1 

0. 0855 

0.009590 

53.09 

0.0778 

0.2071 

35.30 

0.0732 

0.03241 

25.47 

0.0667 

0.04680 

19.97 

0.0617 

0.06214 

16.28 

0.0568 

0.07940 

13.53 

0.0521 

0.09973 

11.43 

0.0472 

0.1233 

9.762 

0.0424 

0.1513 

8.375 

0.0377 

0.1855 

7.179 

0.0329 

0.2288 

6.113 

0.0281 

0.2871 

5.124 

0.0229 

0.3702 

4.083 

0.0169 

0.5211 

2.726 

0.0088 

1.0000 

Since  system  is  closed,  atmospheric  pressure  p.  =  14.696 
psia  has  no  particular  significance  but  serves  as  a  ^convenient, 
if  arbitrary  lower  pressure  limit  for  the  cycle  calculations. 
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PRINCIPAL  RESULTS 


CASE  A 

CASE  B 

CASE 

C  1 

Pe/Al 

KW  SEC 
Kg 

p;/*i 

KW  SEC 
Kg 

nc 

p;/*i 

1400 

0.0028 

3.61 

0.0030 

3.83 

0.0360 

46.14 

1300 

0.0022 

2.86 

0.0041 

5.42 

0.0308 

41.05 

1200 

0.0018 

2.45 

0.0051 

7.06 

0.0260 

36.14 

1100 

0.0013 

1.88 

0.0062 

9.05 

0.0217 

31.40 

1000 

0.0007 

1.09 

0.0072 

10.85 

0.0183 

27.59 

900 

0.0001 

0.15 

0.0080 

12.61 

0.0186 

29.21 

800 

-0.0006 

-0.96 

0.0087 

14.26 

0.0194 

31.82 

700 

0.0092 

15.70 

0.0206 

35.22 

600 

0.0094 

16.81 

0.0220 

39.35 

500 

0.0093 

17.60 

0.0234 

44.17 

400 

0.0086 

17.17 

0.0248 

49.50 

300 

0.0074 

15.69 

0.0257 

54.79 

200 

0.0053 

12.14 

0.0253 

58.27 

100 

0.0018 

4.73 

0.0211 

54.57 

14.696 

-0.0067 

-8.15 

0.0050 

16.65 

Pig.  2A  v  versus  pressure  at  station  1 


Fig.  8A  Net  Electric  Power  versus 
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ELECTRIC  POWER 


2  2  2  2  2 
P.  JeCB  R  p3  A3V3  eCB  R  p3 

e—  ce  ~  J  3  ~  2 

PK3  2P3V3  A3  V3 


Also 


p3  =  Pl  and  V3  =  — 


2  2  2  2 

2  eCB  R  p  eCR  R^p . 

ce  •  “here  - - -  f(hQ  -  hL> 


Reqrouping  factors  - 
2  2 

R  C 
V  =  (— -B,-  ) 


(dimensionless) 


In  metric  units  - 
2  2 

eCB  R\  1 


)  =  4[8 . 854  x  10'12  2^.]  [9.49  x  10 


VTU 


.3  m  °K 
cmb 


x  461 


kg  K 


-  84.73  (j^j) 


In  English  units  - 


eCB2E2 


J_3 

m 


)  =  [84.7 3 [9.480  x  10*4 


x  [3.281%3[0.4536T|j]2 


=  0.5837 


Btu  ft~ 
lbm 


Hence 


v  =  0.5837 


<h0  -  hx) 


(dimensionless ) 
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where 


SUMMARY 

V 


0.5837 


lbm/ft3 

^  (hQ  -  Btu/lbm  J 


ENERGY  BALANCE 


By  a  simple  revision  of  previous  work  we  obtain  the  two  key 
equations . 


-  { cf  +  vy2}  =  0  (2) 


Eq  (2)  can  be  reduced  to  the  form: 

A1{^HHf}  ’  A2{x  +  (~  ~"X)-!}-  A  =  0  (3) 

1  y  x  (y  -  x) 2  J 

where 


2k  — 

2(y  - 

A1 

nE 

(4) 

A2  = 

'4- 

Dy2 

(5) 

A3  " 

+  1  +  cf  +  vy2) 

(6) 

Eq  (3)  can  be  further  reduced  to  the  form: 

B2x2  -  B1x  +  Bq  =  0  (7) 

B2  =  A  +  (3  -  2y)A2  +  A3  (8) 


=  yAx  -  (y2  -  3)A2  +  2yA^ 

b0  =  a2  +  y2a3 

Eq  (7)  may  be  solved  by  the  quadratic  formula, 
of  the  two  roots  is  extraneous  and  may  be  dropped. 

Let  S2  *  B  2  -  4B-B 
1  2  O 


Then  X 


(9) 

(10) 

The  larger 
(11) 

(12) 
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For  the  limiting  case  x  =  1,  Eq  (2)  reduces  to  the  form: 


C2y  -  Cxy  +  CQ  =  0 


where 


-  1  +  v) 


°o  *  +  1  +  °f> 

E 

Let  R2  =  Ci2  -  4C-C 

x  2  0 

The  corresponding  roots  are: 

C1  '  R 
vmin  ~  ^  2C2  ^ 

C1  +  R 

ymax  "  ^  2C~  ^ 

z  .  2 

The  above  roots  are  real  only  if  R  is  not  negative. 
This  will  be  the  case  if  and  only  if : 


c-(-i  -  1)  +  v  (-i  +  1  +  c-)  <  1 
f  nE  nE  f  - 
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